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Abstract 

By generalizing the traditional concept of heat dQ and work dW to also include their time-dependent irreversible components 
diQ and diW allows us to express them in terms of the instantaneous internal temperature T{t) and pressure P(t), whereas 
the conventional form uses the constant values To and Pq of the medium. This results in an extremely useful formulation of 
non-equilibrium thermodynamics so that the first law turns into the Gibbs fundamental relation and the Clausius inequality 
becomes an equality § dQ{t) /T{t) = in all cases, a quite remarkable but unexpected result. We determine the irreversible 
components diQ = diW and discuss how they can be determined to obtain the generalized dW{t) and dQ{t). 
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FIG. 1: Schematic representation of E, E and Eo. We assume 
that E and E are homogeneous and in internal equilibrium, 
but not in equilibrium with each other. The internal fields 
T{t),P(t),- ■ ■ fof E and To,Po,--- of E are not the same 
unless they are in equilibrium with each other. There will be 
viscous dissipation in E when not in equilibrium with E. 



Gislason and Craig [ij recently remarked that the def- 
inition of work in non-equilibrium "...thermodynamics 
processes remains a contentious topic," a rather surpris- 
ing statement, as the field of thermodynamics is an old 
discipline. However, there is some truth to their critique, 
which was motivated by an earlier paper by Bertrand 
a|, who revisited the confusion first noted by Bauman 
3| about different formulation of work dW = PodV or 
dW = PdV in terms of internal (P) and external (Pq) 
pressures [1], see Fig. [U and discussed by many others 
|5l4lOl| since then. Different formulation of work from the 
first law obviously results in different heat. We refer the 
reader to these papers for an interesting history of the 
confusion. Gislason and Craig lal list twenty-six repre- 
sentative textbooks including jlll . Il2| where the pressure- 
volume work and heat are defined so differently that they 
are not equivalent in the presence of irreversibility, and 
there appears to be no consensus about their right for- 
mulation so far [13, p. 181, Vol. 1]. There is obviously no 
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problem for reversible processes. 

The first law of thermodynamics, see for example. Lan- 
dau and Lifshitz [3] or Kondepudi and Prigogine [l5| , 
relates the change dE in the internal energy to heat dQ 
added to and the work dW done by the system 

dE^dQ- dW. (1) 

It is a very general statement and is supposed to be valid 
for all processes. We mostly consider mechanical work, 
but the arguments are valid for all kinds of work; see 
near the end. Zemansky [12I, P-73] defines heat as energy 
exchange "...by virtue of a temperature difference only." 
Unfortunately, this rules out any isothermal (reversible) 
heat exchange and cannot be considered general. Kirk- 
wood and Oppenheim define heat as energy exchange re- 
sulting in "...the temperature increment..." (which rules 
out phase changes requiring latent heat) and later note 
that the work may be converted to heat due to frictional 
dissipation pp. 16, 17] as was first observed by Count 
Rumford in 1798 during the boring of cannon [l5| . There- 
fore, it is natural to account for such viscous dissipation 
in work and heat when dealing with non-equilibrium sys- 
tems. The fact that literature is not very clear on how to 
incorporate viscous dissipation has motivated this work; 
see however [1, [l3| , but the authors do not take the dis- 
cussion far enough to obtain the results derived here. 
Our approach to incorporate viscous dissipation results 
in expressing the first law in terms of internal fields and 
provides an elegant formulation of non-equilibrium ther- 
modynamics in which the Clausius inequality turns into 
an equality in all cases, which is a remarkable result in 
its own right. It has been recently suggested Q that use 
of internal fields is not always consistent with the second 
law. We find no such problem in our approach. 

We consider our system E (see Fig. [1]) surrounded by 
a very large medium S so large that its fields such as its 
temperature Tq, pressure Pq, etc. are not affected by the 
system. They form an isolated system Eq. In the follow- 
ing, all extensive quantities pertaining to E and Eq carry 
an annotation tilde or a suffix 0, respectively, and those 
pertaining to E carry no suffix. Following modern no- 
tation [isl, IJjBI , exchanges with the medium and changes 
within the system carry the suffix e and i, respectively. 
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Traditional Formulation of the First Law To truly ap- 
preciate our contribution, it is useful to consider how the 
first law is traditionally expressed. While the change dE 
is uniquely defined for any infinitesimal process, the value 
of dQ is determined by dW. If dW cannot be uniquely 
determined in irreversible processes because of the am- 
biguity noted above, then dQ will be ill-defined. Tradi- 
tionally, dQ represents the amount of heat exchange dcQ, 
i.e. the heat given to the system from the surrounding 
medium, so that ^^dW is identified with the work ex- 
change —deW = dW = P^dV by the medium to the sys- 
tem, giving deW = PodV since dV — —dV. This is true 
even if the instantaneous pressure P of the system is dif- 
ferent from Pq. As the net heat exchange dcQ + dQ = 0, 
we have deQ = —dQ = —TgdS = TgdcS in terms of the 
entropy change dcS of the system. In other words. 



are given by 14 1 



dE = TdS - PdV, dE = TodS - PodV; 



(3) 



dE = d^Q - d^W = TodeS - PodV 



(2) 



expressed in terms of either exchange quantities or exter- 
nal fields of S. This is Protocol 3 of Bertrand [^]. 

However, the situation is not always clear. Kondepudi 
and Prigogine introduce the work by equating it to PdV, 
where P "...is the pressure at the moving surface," but 
they do not mention whether the form is applicable to all 
processes. Landau and Lifshitz are explicit and state that 
dW = PdV for reversible and irreversible processes 14, 
p. 45]. They require for this the existence of mechanical 
equilibrium (and so do [lll[l3|) within S so that at each 
instant during the process P must be uniform throughout 
the body; its equality with Pq is not required. However, 
they do not discuss dQ when they consider E and E out 
of equilibrium [13, Sect. 20]. If we use dW = PdV for 
the work done by E in terms of P of the system (this 
is similar to Protocol 4 of Bertrand Q), then this will 
alter the heat dQ added to E in Eq. ([Ij for dE must 
be invariant to the choice of internal or external fields. 
To the best of our knowledge, this issue of the actual 
forms of dQ in the two protocols and what is the correct 
form of dW for non-equilibrium systems has not been 
discussed in the literature, even though Kestin [l^ Sect. 
5.12] clearly states that distinguishing heat and work in 
non-equilibrium states is not unambiguous. 

General Consideration and Clausius Relation The 
first law is not useful for any computation unless we can 
ascribe temperatures, pressures, etc. to E and E. This is 
done by taking them to be in internal equilibrium 3, 17 1 
so that their instantaneous entropies are state functions 
S{t) = S{E{t),V{t),---) and S{t) = S{E{t),V{t),- ■ ■) 
of (time-dependent) state variables. They change as re- 
spective state variables E{t), V{t), • • • or E{t), V{t), ■ ■ ■ 
change with time. We first only consider energy and vol- 
ume. The temperatures and pressures are given by ap- 
propriate standard derivatives of the entropies; see Eq. 
([3]) . It now follows that the Gibbs fundamental relations 



quantities without the suffix normally depend on time. 
The validity of Eq. (|3l) requires E and E to be indepen- 
dently homogeneous [14 , 3, 12|- We have dE^ = dVb = 
for Eo. The application of the first law gives 

dE = dQ- dW, dE^dQ- dW, dEo = dQo - dWo = 0. 

We take dQ — d^Q + diQ as a generalization of dQ — d^Q 
in Eq. it denotes the heat added to the system 

either through exchange with its exterior (dcQ) or by 
dissipative internal forces within (diQ). Similarly, dW = 
dcW+ diW = PdV is the work done by the system on its 
exterior (dcW) and by dissipative internal forces (diW). 
A somewhat similar looking approach, but different in 
spirit, is taken in [lil ]. 

To make the our approach computationally useful, we 
need to determine d^Q and diW; see later. With our 
generalization, dQo = diQo and dWo = d-JVo and dQo = 
dWo for Eg. Let us assume that P > Po- The (internal) 
work done by the pressure difference AP = P — Pq > 
is dWo = APdV > , since dV > 0. This resuhs in 
raising the kinetic energy dKg of the center-of-mass of the 
surface separating E and E and overcoming work dWfr 
done by all sorts of viscous or frictional drag. Thus, 

dWo = diWo = dKs + dWh. 

Because of the stochasticity associated with any statis- 
tical system, both energies on the right side dissipate 
among the particles so as to increase the entropy and 
appear in the form of heat (= c?iQo > 0) within the iso- 
lated system. This is the generalized heat dQo > 0, since 
dcQo = for Eq. It is normal to associate all irreversible 
components in Eq with E {diQ = diQo,diW = diWo) 
and not with the (extensively large) medium E, which 
always remain in internal equilibrium [itI ] ; compare with 
the situation of finite surroundings considered by Bizarro 
Q. Thus, when there are irreversible processes going on, 
it is natural to generalize heat from dcQ in Eq. ([2]) to 
include the irreversible heat diQ and identify dQ as the 
heat added to the system. Similarly, we need to general- 
ize work from dgW = PodV to dW = PdV and identify 
it as work done by the system. Remarkably, we find that 

diQ = diW. (4) 
The generalization docs not change dE and we have 

dE = d^Q - d^W = dQ-dW ^dQ- PdV. (5) 
However, the most important result is that [l3] 

dQ = TdS, (6) 

even when the system is not in equilibrium with the 
medium. This is easily seen by comparing Eq. ([5]) with 
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the Gibbs fundamental relation in Eq. ([3]). This makes 
dQ/T an exact differential dS so that 



j dQ/T : 



for any cyclic process. (7) 



This leads to 



,Q/T = -§d,Q/T > 0, (8) 



which is consistent with the traditional Clausius inequal- 
ity §dcQlT < 0. Thus, the Clausius equality in Eq. 
should not be interpreted as the absence of irreversibility, 
as is clear from the above inequality. It is only because 
of the use of the generalized heat dQ that this inequal- 
ity has become an equality. The second integral with its 
sign in Eq. ^ can be thought of as the uncompensated 
transformation N [15]. Our formulation has allowed us 
to identify it with the first integral 



N = §diQ/T - fdiWo/T, 



(9) 



which provides a way of computing N using equations of 
state of E and S. From diS = dS — dcS, we find that 

§diS = -§d,Q/To > 0, 

which is the second law for a cyclic process and represents 
the irreversible entropy generated in a cycle; recall that 
dcQ is the traditional dQ in Eq. 

General Form of Work Let us follow the conse- 
quences of this particular extension a bit further and 
prove that only dW = PdV is consistent with the^econd 
law. Since dQo = dQ + dQ = dWo = dW + dW > 0, 
we find that the net heat dQo in the medium and the 
system need not vanish, in contrast to what happens in 
equilibrium: 



dQ + dQ^ diWo > 



(10) 



because of viscous dissipation diQa = diWo, and is re- 
lated to the irreversible entropy generation within Sq 
or E, as we will see later. For frictional forces only 
(P = Pq) which are included in our approach [l9| , the 
above relation reduces to Eq. (7) of Bizarro [8|. From 
dW + dW = PdV + PodV valid for any arbitrary dV, it 
follows that 



dW = PdV, dW = P(,dV, 



(11) 



which proves the above assertion. The alternate choice 
dW = PodV and dW = PdV employed in the tradi- 
tional formulation, see Eq. ([2|), will violate the second 
law as the choice will lead to negative heat dissipation 
dQo = c?Wo = {Pq — P)dV < 0, a physical impossibility. 
Thus, we must write the first law for the system and the 
medium, respectively, as 



The above demonstration establishes that the work is 
given by Eq. (jlip in terms of the instantaneous pressure 
in all cases contrary to the traditional wisdom, see Eq, 
([2]), that it is given by PodV. The generalized formula- 
tion brings out a symmetry in the form of dW and dW 
under the interchange system<^=^>medium: the work for 
any body is always uses its own internal pressure. This 
symmetry is absent in the traditional formalism in which 
both terms are given by the external pressure Pq- 

Benefits of the Generalized Formulation The new for- 
mulation has many other desirable properties. For ex- 
ample, the entire thermodynamics and stability consid- 
erations for any body can be expressed in terms of the 
variables associated with the body alone at each instant. 
We can use the equations of state of the body alone for 
thermodynamic computation. The first law is no longer 
different from the Gibbs fundamental relation so we only 
deal with equalities and not inequalities. The general- 
ized heat dQ and work dW differ from dcQ and dcW, 
respectively, by the same contribution diW = diWo; the 
latter is the term dT, in Eq. (19) of Eu \2C^, which can be 
absorbed in his work term (= —d^W in our notation) to 
reduce this equation to the above form for dE. The de- 
termination of diWo is straight forward by measuring P 
and Pq. Adding this to dcQ then allows us to determine 
dQ. Consider, as an example, the case when the medium 
pressure Pq = such as when a gas expands in vacuum. 
In this case, dgW = 0; dW is the irreversible work, which 
cannot be less than zero, as it is given by the irreversible 
entropy generation as shown below. 

Irreversible entropy generation in an isothermal pro- 
cess We now give an alternative derivation of Eq. (jlip 
by computing the irreversible entropy gain so that the 
role of the latter can be clearly seen in an irreversible 
process. A brief discussion of this issue can be found 
in Kondepudi and Prigogine [11, pp. 94-95] where they 
discuss irreversible expansion of a gas; see also [l3]- We 
will assume that the temperature of the system and the 
medium is constant during the process at Tq. The pres- 
sure difference gives rise to dissipation, which results in 
an irreversible heat generation diQ — TodiS in addition 
to the heat deQ = —dQ = —TgdS — TodeS received 
from the mediurn; here, diS — diS''-^-' is the irreversible 
entropy gain 15|, [iTf 



diS'(^) = iP- Po)dV/To > 0. 

We see that dQ = To{dcS + diS) — TodS in accordance 
with Eq. ^ so that by including diQ in dQ relates dQ to 
the entropy change dS. The generalized first law for the 
medium in Eq. (fT2)) reads dE = TodS — PodV, in which 
we can replace dV by —dV and dE by —dE to obtain 
dE = —dQ — PodV. We compare this with the right side 
of Eq. H]) for the system and use the fact that 



dE^dQ- PdV, dE^dQ- PodV. 



(12) 



dQ + dQ — diQo — d^Q 
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to determine the work dW. Simple algebra immediately 
leads to dW = PdV, as claimed above. We also note that 
we do not need Jarzynski result to obtain diS^^\ 

General Irreversible Process We now consider the 
same process as above, except that the temperature T 
of the system is also allowed to be different from that 
of the medium, which we take to be the constant Tq. 
The general relation in Eq. (|4]) is still valid. However, 
diS consists of diS^^'> and diS'(Q) or diS'^^); the former is 
given above and the latter contribution is given by 

diS'(Q) = (l/T-l/To)rfQ = -{T-To)dS/To = diS'^^' > 0. 

We immediately find that diS'*^^-' and diS^^^ are formally 
of the same form; see 17[ . To be convinced of the form of 
(iiS'(Q) = diS'(^\ we start with the identity d^Q-PadV = 
TdS — PdV , which can be transformed into TodiS ~ 
(Tq - T)dS + {P - Po)dV We can now identify 

the above two components of diS. To obtain diS''-*^', we 
use Eq. (|6]) to replace dS. Both irreversible components 
are non-negative as expected from the second law. We 
thus see that expressing the first law in terms of internal 
fields is consistent with irreversible entropy production, 
contrary to a recent claim Q . 

Both irreversible entropy gains can be easily deter- 
mined. The determination requires measuring the two 
temperatures and pressures, while the determination of 
dQ has already been explained. Thus, the generalized 
formulation creates no additional problem in experimen- 
tally determining various quantities. 

Inclusion of other state variables Let us now extend 
the discussion to include other extensive quantities such 
as the flow of matter, the electric interactions, chemi- 
cal reactions, etc. For specificity, we focus on chemical 
reactions, which we assume to be described by a single 
extent of reaction ^. The corresponding affinity for the 
system is given by A, while that for the medium is given 
hy Aq — 0. We assume another observable X such as the 
number of solvent in a binary mixture. The correspond- 
ing chemical potential is /x for the system and /io for the 
medium. The work is now dW = PdV - fidX + Ad£, [H | . 
The Gibbs fundamental relation for E is given by dE = 
TdS-PdV+ndX-Ad^, while the first law for it takes the 
form dE — dQ — PdV + fidX — AdS,. Comparing the two 
we find that Eq. (O still holds so that the validity of Eq. 
([7]) is not affected by the presence of other variables in the 
first law Rewriting dQ = dE + PdV - ^idX + Ad£, as 
dQ = dE+PodVpjfidX+{P-Po)dV-{iJ,-no)dX+Ad(, 
we can identify [19| deQ with the first three terms in 
which deW = PodV — nodX for work of all types [2l[. 
Similarly, diQ = d^W = dW^ represents the last three 
terms (P - Po)dV - (a* - tJ-o)dX + AdS,- see Eq. (g]). We 
again obtain the general results in Eq. ([5]), as expected. 
According to the second law, each contribution in d\W 
is non-negative. The extension to arbitrary number of 
observables and internal variables is trivial 



In summary, we have shown that the definition of the 
generalized heat dQ, which includes its irreversible com- 
ponent diQ, follows uniquely from the unique choice of 
dW resulting from the second law; see Eq. (fTT|) . It in- 
cludes its irreversible component diW = d^Q; see Eq. 
(PI). Another remarkable consequence is that in terms 
of generalized dQ, the Clausius equality ([7]) is always 
maintained, in contrast to the inequality ^ in the tra- 
ditional approach. Our generalized formulation brings 
about a symmetry between the system and its surround- 
ing medium, see Eq. pip , which is absent in the tradi- 
tional approach using external fields. The first law be- 
comes identical with the Gibbs fundamental relation so 
that we only deal with equalities that are easier to deal 
with than the inequalities in the traditional approach. 
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